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Abstract. We consider the linear Wigner-Fokker-Planck equation subject to 
confining potentials which arc smooth perturbations of the harmonic oscillator 
potential. For a certain class of perturbations we prove that the equation ad- 
mits a unique stationary solution in a weighted Sobolev space. A key ingredient 
of the proof is a new result on the existence of spectral gaps for Fokker-Planck 
type operators in certain weighted L'^— spaces. In addition we show that the 
steady state corresponds to a positive density matrix operator with unit trace 
and that the solutions of the time-dependent problem converge towards the 
steady state with an exponential rate. 



1. Introduction 

This work is devoted to the study of the Wigner-Fokker-Planck equation (WFP), 
considered in the following dimcnsionlcss form (where all physical constants are 
normalized to one for simplicity): 

dtw + £_ ■ \7xw + e[V]w = A^w + 2 div^ {£_w) + A^w, 

where x,£, G M'', for d> 1, and t G K+. Here, w{t,x,£,) is the (real valued) Wigner 
transform j37j of a quantum mechanical density matrix p{t,x,y), as defined by 

^ f f V 

(2^)° 

Recall that, for any time t G ]R-|_. a quantum mechanical (mixed) state is given by a 
positive, self-adjoint trace class operator pit) G 3/^^ . Here we denote by ^{L'^{M.'^)) 
the set of bounded operators on L^(M'^) and by 

^1 := {p e m{L'^{M.'^)) : tr \p\ < oo} , 

the corresponding set of trace-class operators. We consequently write p e £^1^ C 
,3^, if in addition p > (in the sense of non- negative operators). Since C ^2, 
the space of Hilbert- Schmidt operators^ i.e. 

:= {p e ^(L2(M'i)) : tr(p» < 00} , 



(1.1) 



(1.2) .(,.,0:^---,/ p(t,. + |,.-|)e-^-"d. 
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we can identify the operator p{t) with its corresponding integral kernel p{t, ■,■) G 
the so-called density matrix. Consequently, p{t) acts on any given function 

(p e yia 

{p{t)ip){x) = j p{t,x,y)ip{y)dy. 

Using the Wigner transformation (jl.2|) . which by definition yields a real-valued 
function w{t,-,-) € L^(R^'^), one obtains a phase-space description of quantum 
mechanics, reminiscent of classical statistical mechanics, with x G R"* being the 
position and ^ G M'' the momentum. However, in contrast to classical phase space 
distributions, w{t,x,^) in general also takes negative values. 

Equation (jl.ip governs the time evolution of w{t,x,£,) in the framework of so- 
called open quantum systems, which model both the Hamiltonian evolution of a 
quantum system and its interaction with an environment (see |13| . e.g.). Here, 
we specifically describe these interactions by the Fokkcr-Planck (FP) type diffusion 
operator on the r.h.s. of (jl.ip . For notational simplicity we use here only normalized 
constants in the quantum FP operator. However, all of the subsequent analysis 
also applies to the general WFP model presented in [34] (cf. Remark 12.41 below). 
Potential forces acting on w{t, •, •) are taken into account by the pseudo-differential 
operator 

(1.3) ie[V]f)ix, := --i-^ // 6V{x, v) fix, O e'^<^-^"> d^' dr,, 
where the symbol SV is given by 

(1.4) SV{x,7j)^v[x+^)-v[x~^), 

and 1^ is a given real valued function. The WFP equation is a kinetic model 
for quantum mechanical charge-transport, including diffusive effects, as needed, 
e.g., in the description of quantum Brownian motion [TS], quantum optics [T5], 
and semiconductor device simulations |16j . It can be considered as a quantum 
mechanical generalization of the usual kinetic Fokker-Planck equation (or Kramer's 
equation), to which it is known to converge in the classical limit ?i — > 0, after an 
appropriate rescaling of the appearing physical parameters |10| . The WFP equation 
has been partly derived in (llj as a rigorous scaling limit for a system of particles 
interacting with a heat bath of phonons. Additional "derivations" (based on formal 
arguments from physics) can also can be found in [HJ [TSl 133 |3S • 

In recent years, mathematical studies of WFP type equations mainly focused on 
the Cauchy problem (with or without self-consistent Poisson-coupling) , see [2J |31 
131 [3 m [12]. In these works, the task of establishing a rigorous definition for the 
particle density n(t, x) has led to various functional analytical settings. To this 
end, it is important to note that the dynamics induced by (|l.ip maps =^+(L^(R'^)) 
into itself, since the so-called Lindblad condition is fulfilled (see again Remark 12.41 
below). For more details on this we refer to [7] [3] and the references given therein. 
In the present work we shall be mainly interested in the asymptotic behavior as 
t — >■ -l-oo of solutions to (jl.ip . To this end, we first need to study the stationary 
problem corresponding to (|l.ip . Let us remark, that stationary equations for open 
quantum systems, based on the Wigner formalism, seem to be rather difficult to 
treat as only very few results exist (in spite of significant efforts, cf. |B] where the 
stationary, inflow-problem for the linear Wigner equation in d = 1 was analyzed). 
In fact the only result for the WFP equation is given in [34], where the existence 
of a unique steady state for a quadratic potential V{x) oc |xp has been proved. 
However, the cited work is based on several explicit calculations, which can not be 
applied in the case of a more general potential V{x). 



STATIONARY STATES AND LARGE TIME BEHAVIOR OF WFP 



3 



The goal of the present paper is twofold: Fhst, we aim to estabhsh the existence 
of a normahzed steady state Wooixj^,) for in the case of confining potentials 
V{x), which are given by a suitable class of perturbations of quadratic potentials 
(thus, V{x) can be considered as a perturbed harmonic oscillator potential). The 
second goal is to study the long-time behavior of p.ip . We shall prove exponential 
convergence of the time-dependent solution towards Wca as t — > +00. 

In a subsequent step, we shall also prove that the stationary Wigner function Woo 
corresponds to a density matrix operator Poo G ■ Remarkably, this proof exploits 
the positivity preservation of the time- dependent problem (using results from [9]), 
via a stability property of the steady states. 

To establish the existence of a (unique) steady state Woo, the basic idea is to 
prove the existence of a spectral gap for the unperturbed Wigner-Fokker-Planck 
operator with quadratic potential. This implies invertibility of the (unperturbed) 
WFP-operator on the orthogonal of its kernel. Assuming that the perturbation 
potential is sufficiently small with respect to this spectral gap, we can set up a 
fixed point iteration to obtain the existence of Woo- The key difficulty in doing 
so, is the choice of a suitable functional setting: On the one hand a Gaussian 
weighted L^-space seems to be a natural candidate, since it ensures dissipativity of 
the unperturbed WFP-operator (see Section[3]). Indeed, this space is classical in the 
study of the long-time behavior of the classical (kinetic) Fokker-Planck equation, see 
[26| . However, it does not allow for feasible perturbations through 6[Vo]. In fact, 
even for smooth and compactly supported perturbation potentials V07 the operator 
8[Vo] would be unbounded in such an L^-space (due to the non- locality of 8[Vo], see 
Remark 1 5. 2p . We therefore have to enlarge the functional space and to show that 
the unperturbed WFP-operator then still has a (now smaller) spectral gap. This is 
a key step in our approach. It is a result from spectral and semigroup theory (cf. 
Proposition l4.8p which is related to a more general mathematical theory of spectral 
gap estimates for kinetic equations, developed in parallel in [M] (see also [IS]). We 
also remark that for V{x) = the WFP equation corresponds to a differential 
operator with quadratic symbol [34j and thus our approach is closely related to 
recent results for hypo-elliptic and sub-elliptic operators given in [17l [26l [31] . 

Comparing our methods to closely related results in the quantum mechanical 
literature, we first cite |20| . where several criteria for the existence of stationary 
density matrices for quantum dynamical semigroups (in Lindblad form) were ob- 
tained by means of compactness methods. In [S] the applicability of this general 
approach to the WFP equation was established. In [22l [21] sufficient conditions 
(based on commutator relations for the Lindblad operators) for the large-time con- 
vergence of open quantum systems were derived. However, these techniques do not 
provide a rate of convergence towards the steady states. In comparison to that, 
the novelty of the present work consists in establishing steady states in a kinetic 
framework and in proving exponential convergence rates. However, the optimality 
of such rates for the WFP equation remains an open problem. In this context one 
should also mention the recent work |25] . in which explicit estimates on the norm 
of a semigroup in terms of bounds on the resolvent of its generator are obtained, 
very much along the same lines as in present paper and in |24j . 

The paper is organized as follows: In Section[2]we present the basic mathematical 
setting (in particular the class of potentials covered in our approach) and state our 
two main theorems. In Section [3] we collect some known results for the case of a 
purely quadratic potential and we introduce the Gaussian weighted L^-space for 
this unperturbed WFP operator. This basic setting is then generalized in Section[?l 
which contains the core of our (enlarged) functional framework: We shall prove 
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new spectral gap estimates for the WFP operator with a harmonic potential in 
spaces with only polynomial weights. In Section [S] we prove the bomidedness of the 
operator 0[Vb] in these spaces. Finally, Section |6] concludes the proof of our main 
result by combining the previously established elements. Appendix \K\ includes the 
rather technical proof of a preliminary step which guarantees the applicability of 
the spectral method developed in [53]. 

2. Setting of the problem and main results 

2.1. Basic definitions. In this work wc shall use the following convention for the 
Fourier transform of a function /^{x): 

^{k) := [ e-'^^-^dx. 

From now on we shall assume that the (real valued, time-independent) potential 
V, appearing in is of the form 

(2.1) Vix) = ^\x\^ + XVoix), 

with Vq e C°°(]R'';]R) and A G M some given perturbation parameter. In other 
words we consider a smooth perturbation Vq of the harmonic oscillator potential. 
The precise assumption on Vq is listed in (|2.9p . An easy calculation shows that for 
such a V the stationary equation, corresponding to (jl.ip . can be written as 

(2.2) Lw = Xe[Vo]w, 
where L is the linear operator 

(2.3) Lu; • W^w + x ■ V^w + A^w + 2 div^ {£^w) + A^w. 
Remark 2.1. When considering the slightly more general class of potentials 

V{x) = ^\x\^ + a-x + XVo{x), XeR,a<EM.'', 

we would find, instead of (j2.3p . the following operator: L^w Lw + a ■ V^w. 
Thus, by the change of variables x t-^ x + a wc are back to (j2.3p . 

The basic idea for establishing the existence of (stationary) solutions to (|2.2p is 
the use of a fixed point iteration. However, L has a non-trivial kernel. Indeed it 
has been proved in [33] that, in the case A = 0, there exists a unique stationary 
solution /I e iS(M^'^), satisfying 

(2.4) Lfi = 
and the normalization condition 

(2.5) // ti{x,Odxd^^l. 

Explicitly, /i can be written as 

(2.6) M-ce-^("'«), 
where the function A is given by 

(2.7) A{x,0--^\{\^\' + 2x-( + 3\e), 

and the constant c > is chosen such that (|2.5p holds. Note that for any p G 
such that w £ L^(R^'*) the following formal identity 

ti p = / p{x,x)dx= // w{x,£_) dx d^,, 
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can be rigorously justified by a limiting procedure in see [T]. Since trp is 
proportional to the total mass of the quantum system, we can interpret condition 
(|2.5|) as a mass normalization. 

In the following, we shall denote by ct > the biggest constant such that 

(2.8) Hcss^ - cr I > 0, for all {x, ^) G M^'', 

in the sense of positive definite matrices, where I denotes the identity matrix on 
M^''. In the analysis of the classical FP equation, condition (|2.8p is referred to as 
the Bakry-Emery criterion [8] . In our case one easily computes 

cr = 1 - 1/V2. 

In a Gaussian weighted L^-space, a will be the spectral gap of the unperturbed 
WFP-operator and hence the decay rate towards the corresponding steady fi (cf. 
(EH), (|XTU|) below). 

The functional setting of our problem will be based on the following weighted 
Hilbert spaces. While the stationary and transient Wigner function are real valued, 
we need to consider function spaces over C, for the upcoming spectral analysis. 

Definition 2.2. For any m G N, we define Hm := ^^(M^'', t/^Ma; d^), where the 
weight is 

We equip 'Hm with the inner product 

(/'•9>«„, = // — dxde 
Clearly, we have that Hm+i C Hm, for all m G N. 

2.2. Main results. With these definitions at hand, we can now state the main 
theorems of our work. Note that for the sake of transparency we did not try to 
optimize the appearing constants. 

Theorem 1. Let m > Kd he some fixed integer, where K = K{A) G (1, 144] is 
a constant depending only on A{x,(^) (defined in Lemma \4.2^ . Assume that the 
perturbation potential Vq satisfies 

(2.9) r,„ := Cm. max ||5;JVo|1l=-(R'') < +00, 

\j\<m 

where Cm > depends only on m and d, as seen in the proof of Proposition I5.il 
Next we fx some 7™ G (0,7m), where j„i > is given in (j4.7D . Furthermore, let 
the perturbation parameter X satisfy 

(2.10) 1^1 <-F^' 

where 5m = <5m(7m) > 1 is defined in (j4.15|) . Then it holds: 

(i) The stationary Wigner- Fokker- Planck equation (|2.2p admits a unique weak 
solution Woo G Hm H H^{M.'^'^), satisfying //]g2d li'oo dx d^ ~ 1. Moreover, 
Woo is real valued and satisfies iUqo S -^locC^^'')- 

(ii) Equation (|l.ip admits a unique mild solution w G C([0, c»), "Hm)- In addi- 
tion, for any such mild solution w{t) with initial data wq G Hm satsifying 
JJs2d wq dxd^ = 1, we have 

II W{t) - Woo ||^^,„ < Sme'""-* II Wo - Woo 11^^ , Vi > 0, 

with an exponential decay rate 

Km := 7m - |A|(5m.rm. > 0. 
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(iii) Concerning the continuity of Woo ~ Woo{^) w.r.t. X, we have 

11 II ^ \X\OmTm ,1 II 

II Woo " M ||^^„ < ~ r ll^ll^-- 

Remark 2.3. In this result, the constant (t,„ :~ ^m/6m > 0, roughly speaking, 
plays the same role for L on Hm as cr > docs in the case of H (where H is defined in 
Definition [5TT]), where it is nothing but the size of the spectral gap, see Proposition 
13.51 For L on "H, cr also gives the exponential decay rate in the unperturbed case 
A = 0. For L on Hm the situation is more complicated. Here, assertion (ii) yields 
an exponential decay of the unperturbed semi-group with rate Km = 7m G (0,7m) 
and 7„i ^ am (but possibly equal to cr, as can be seen from (|4.15[) '). In addition, 
one should note that 5m > 1 may blow-up as 7„i 7„i, cf. estimate ()4.1ip . 

Theorem [1] is formulated in the Wigner picture of quantum mechanics. We shall 
now turn our attention to the corresponding density matrix operators p(t). This 
is important since it is a priori not clear that Wao is physically meaningful - in 
the sense of being the Wigner transform of a positive trace class operator. To this 
end we denote by poo the Hilbert-Schmidt operator corresponding to the kernel 
Poa{x,y), which is obtained from Woo{x,S,) by the inverse Wigner transform, i.e. 

Analogously we denote by po the Hilbert-Schmidt operator corresponding to the 
initial Wigner function G T-Lm- 

We remark that the existence of a unique mild solution of equation p.ip on T-Lm 
will be a byproduct of our analysis. 

Theorem 2. Let m > Kd he some fixed integer. Let Vo, A, and wq satisfy the 
same assumptions as in Theorem]^ Then we have: 

(i) The steady state poo *s o, positive trace-class operator on L^(M'*), satisfying 
trpoo = 1- 

(ii) Let p G C([0,oo), ^) be the unique density matrix trajectory corresponding 
to the mild solution of (jl.ip . Then, the steady state poo is exponentially 
stable, in the sense that 

II p{t) - Poo 11^, < (27r)^5™e-""*|| Wo - Woo 1^^ , > 0. 

(iii) // the initial state wq G Tim corresponds to a density matrix po G ( and 
hence wq is real valued, tr pQ = JJ wq dxd£, = 1 /, then we also have 

lim II p{t) - PooW^r, = 0. 

Note that, in the presented framework, we do not obtain exponential convergence 
towards the steady state in the i?i-norm but only in the sense of Hilbert-Schmidt 
operators. This is due to the weak compactness methods involved in the proof of 
Griimm's theorem (cf. the proof of Th. [2]in 

Remark 2.4. Consider now the following, more general quantum Fokker-Planck 
type operator replacing the r.h.s. of 

Qw :— DppA^w + 2 Dpq divx (V^w) + 2Di div^ (^w) -f- Dqq A^w. 

It is straightforward to extend our results to this case as long as the Lindblad 
condition holds, i.e. 

(2.11) Dpp > 0, DppD^^ - (dI^ + ^) > 0. 
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The modified quadratic function A{x,^) is given in |34j . The Lindblad condition 
()2.1ip implies tliat discarding in (jl.ip tlic chffusion in x, and hence reducing the 
r.h.s. to the classical Fokker-Planck operator Qc\w := A^w + 2 div^(^u;), would not 
describe a "correct" open quantum system. Nevertheless, this is a frequently used 
model in applications [35]: yielding reasonable results in numerical simulations. 

3. Basic properties of the unperturbed operator L 

3.1. Functional framework. It has been shown in [31] that the operator L, de- 
fined in (|2.3p . can be rewritten in the following form 

(3.1) Lw = div (Vw + w(VA + i^)) , 
with 

(3.2) div(Fe-^) = -div(F/i) = 0. 

c 

Here and in the sequel, all differential operators act with respect to both x and ^ 
(if not indicated otherwise). In p.ip . the function A is defined by (|2.7p and 

The reason to do so is that p.ip belongs to a class of non-symmetric Fokker-Planck 
operators considered in [8] . From this point of view, a natural functional space to 
study the unperturbed operator L is given by the following definition. 

Definition 3.1. Let H := ^^(IR^'^,^-! dxd^), equipped with the inner product 

{f^9)u- [[ 

We can now decompose L into its symmetric and anti-symmetric part in i.e. 

(3.4) L = L'' + L^\ 
where 

(3.5) ~ div (Vw + wVA) , L''^w:= div (Fw). 
It has been shown in [53], that the following property holds: 

(3.6) L^'fi = 0, L^'a^ = 0. 

where n is the stationary state defined in (|2.6p . Next we shall properly define 
the operator L. To this end we first consider L\„^, which is closable (w.r.t. the 

■H-norm) since it is dissipative: 

Lemma 3.2. i|poe «s dissipative, i.e. it satisfies }ic{Lw,w)^ < 0, for all w G 
C^(R2'*). 

Proof. Using VA = — /i^^V/i we have, on the one hand 

{Vw,w)n,^ II — div (Vw -f wV A) dx dC = // — div f ( — ) ) d.T d^ 

J JKi-i M J Jv?i A* V \l^J J 

1 

dx di < 0. 

On the other hand, it follows from p.2[) that 

w 

wdivF = F- V^, 

M 

and thus 

div(i^^.) = -f.F .(^v^,-^)^^^F■v(^ 
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An easy calculation then shows 



Rc{L'^w,w)-^=Rc 1 1 -dw{Fw)dxd^^RcJJ - F ■ V ] fidx 



1 



/ 


W 







2 , , div(F/i) da; = 0, 

by p.2[) . To sum up we have shown that Re {Lw, w)^ < holds. □ 



The operator L = L\„^ is now closed, densely defined on "H and dissipative. 
Moreover one easily sees that L* = — L'^^. The main goal of this section is to 
prove that L admits a spectral gap and is invertible on the orthogonal complement 
of its kernel. For the first property, we start showing that L is the generator of a 
Cp-semigroup of contractions on H. For this we recall the following result from [9]. 

Lemma 3.3. Let the operator P — p2 (a;, ^, V^, V^), where p2 is a second order 
polynomial, be defined on the domain &{P) = C^{M.'^'^). Then P is closable and 



is the maximum extension of P in L'^{M.'^'^) 



Several variants of such a result (on different functional spaces) can be found in 
[3] . We can use this result now in order to prove that L is the generator of a 
Co-semigroup. 

Lemma 3.4. L generates a Co-semigroup of contractions on %. 
Proof. Defining v :~ w/ ^JJl transforms the evolution problem 

dtw = Lw, ^\i-o = ""^0 G ^ 

into its analog on L'^{^^'^). The new unknown f(i,x,^) then satisfies the following 
equation 

dtV = Hv, wo/y/JI, 
where H is (formally) given by 

Hv = Av + F - Vv + Uv, 



and the new "potential" U ~ U{x,^) reads 

U 

Defining H on ^(H) = C^{R'^'^), we have 



U = ^AA-^\VA\'. 



Lw = \/TiH 



w \ 



and thus the dissipativity of L on Ti directly carries over to H = Hi on L (I 
Next, we consider H*\ , defined via {Hf,g)r2 — {f,H*g)r2, for f,g € C^(I 
Due to the definitions (12. 7p and p. 31) , the operator H* I „^ is exactly of the form 

needed in order to apply Lemma 13.31 Thus, H* = H*\^^ is also dissipative (on all 

" __ 

of its domain). Hence, the Lumer-Phillips Theorem (see [30]! Section 1.4) implies 
that H is the generator of a Co-semigroup on L'^iM?'^), denoted by e^*. 

Reversing the transformation w ^ v then implies that L is the generator of the 
Co-semigroup Ut onH, given by 



This finishes the proof. □ 
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3.2. Semigroup properties on H. The above lemma shows that the unperturbed 
WFP equation 

dtw = Lw, If Ij^q ~ wq gT-L, 

has, for all wq G H, a. unique mild solution w g C([0, oo), 'H), where w{t,x,^) = 
UtWo{x,^), with Ut defined above. Obviously we also have Utfi = fi, by (|2.4p . 
Moreover, in [53] the Green's function of Ut was computed explicitly. It shows that 
Ut conserves mass, i.e. 

// w{t, X, ^)dxdi ^ // wo{x,Odxd^, Vi > 0. 

JJK2d J jR^d 

Next, we define 

(3.7) ■.= {weH■.wJ.^i}cH, 

which is a closed subset of H. Note that w J- fi simply means that 

{w,^)^= Ij w{x,£^)dxd^ = Q. 

Hence, we have for w g C^(R^''), using (jSTS)) : 

{L'''w,fi)^= [[ L'''wix,C)dxd^ ^0. 



Thus, L'^" : n-^ n ^(L^") ^ n-^. Moreover, : V.-^ n &{L^) , since 

is spanned by the eigenfunctions of (except of fi). To sum up, the operators 
and L^^ arc simultaneously reducible on the two subspaccs T-L = span[/i] ® Ti,^ . 

We also have that Ut maps into itself, since for wq G the conservation 
of mass implies 

(3.8) (UtWQ,^i)-^= w{t, X, ^)dxd^ ^ wo(x, C) dx d^ = 0, > 0. 

Lemma allows us to prove that L has a spectral gap in H, in the sense that 

(3.9) cr(L) \ {0} C {z e C : Rez < -cr}. 
Proposition 3.5. It holds 



where a > is defined in 



Proof. Condition (|2.8p implies that has a spectral gap of size ct > (cf. §3.2 in 
[S], e.g.). Moreover, [51 Theorem 2.19] also yields exponential decay (with the same 
rate) for the non-symmetric WFP equation: 

(3.10) II Ut{wo - m)II« < c^"*|| wa - fi \\^. 

Here, wq G H has to satisfy JJ^2d wq dxd^ = JJ^2d M dxd^ = 1. By the discussion 
above, we know that L\^^ is the generator of Ut\^_^- Hence, p.lOp implies 

(3.11) ||(L-z)-^||,^(^^) < VzeC, Rez>-a, 

which proves the assertion for z = 0. □ 

As a final preparatory step in this section, we shall prove more detailed coercivity 
properties of L within H-^. We shall denote := {w € H : Vw S H}, and 'H~^ 
will denote its dual. 
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Lemma 3.6. In Ti.^ the operator L satisfie 



(3.12) 



Rc {Lw, w)^ > j 



Similarly, there exists a constant Q < a < a , such that 

(3.13) -Re {Lw,w)^>a\\w\\\i, WweH^nn^. 

Proof. Wc shall use the weighted Poincare inequality (see [8]): For any function 
/ e L^{M:^'',i.idxd^), such that JJ^^^ f^J.dxdC = 0, it holds: 



(3.14) 



Estimate p.l2p then readily follows by setting / = w/^: 

2 



Re {Lw, w)u 



R2d 



dx d^ < —a 



R2d ^ 



dxdS,. 



In order to prove assertion (|3.13p . we note that 



-2d- 



— div 



taking into account that A(log/i) = — 2(i. Next, let < a < 1 (to be chosen later), 
and write 

2 



Re {Lw, w)-H = — 



V , 



dxd£. -{1-a) 



(1-a) 



dx dC + 2da 

2 



W 



dardC 



dxdC 



dxd^. 



Inequality (|3TT4)) for / 

-(1-a) 
Therefore 

Re {Lw, w)u < —a 



ui//^ then implies: 

2 



dxdC < -f7(l - a) 



dxd£,. 



dxdS, + {2da-(j{l-a)) jj ^— ^ dxd^. 

[{2d ^ 



The choice a = a/{a + 2d+l) yields assertion p.l3p . 



□ 



In the next section we shall study the operator L in the larger functional spaces 
V-m (see Definition [22]). This is necessary since the perturbation operator 0[Vb] is 
unbounded in %, even for Vb G C^iW^), cf. Remark[ 



4. Study of the unperturbed problem in 



In this section, we adapt the general procedure outlined in [24l |29] to the specific 
model at hand. One of the main differences to the models studied in is the fact 
that the WFP operator includes a diffusion in x. Nevertheless, we shall follow the 
main ideas of [23]. In a first step, this requires us to gain sufficient control on the 
action of Ut on T-Lm- After that, we establish a new decomposition of L (not to be 
confused with the decomposition L — L^ + L^^ used above) in order to lift resolvent 
estimates onto the enlarged space Hm 3 "H. Together with the Gearhart-Priiss 
Theorem (cf. Theorem V.1.11 in [H]), these estimates will finally allow us to infer 
exponential decay of Ut on Jim. 



STATIONARY STATES AND LARGE TIME BEHAVIOR OF WFP 



11 



4.1. Mathematical preliminaries. In p.Sp we decomposed the unperturbed evo- 
lution operator as i = L'* + L^^. As a first, from basic property of the spaces Tim 
(see Definition 12. 2[) , we note that L'^ is still anti-symmetric in Hm , m € N. 

Lemma 4.1. It holds 

(4.1) Re{L''''w,w)-H,„ =0, Vm G N. 

Proof. A straightforward calculation yields div(F) = 0. Hence 

Rc{L'''w,w)^ ^Rc [[ —div{Fw)dxd^= ^ [[ F ■ \/ {\w\^) i^-^ dx 



\w\^ m A"^-^ F - V Adxd£„ 



after integrating by parts and using v^-^ = 1 + A™(x, ^). Now, it is easily seen from 
(1331) that F • VA = 0, which implies gJ]). □ 

The proof shows that, in the definition of "Hm , it is important to choose the weight 
h^'tyi, as a (smooth) function oi A = — log ^ . Otherwise the fundamental property 
(|4.ip would no longer be true. Also note that in contrast to L^^, the operator is 
not symmetric in 'Hm- Before studying further properties of L in Tim we state the 
following technical lemma. In order to keep the presentation simple, we shall not 
attempt to give the optimal constants. 

Lemma 4.2. Let A = — log^, as given in (|2.7p . Then the following properties 
hold: 

(a) There exists a constant ai > 0, such that for all m Cz N it holds: 

ai (1 + A™) < ^"-1 \VAf , for all \x\'^ + \^\^ > 12. 

(b) Choosing K :— it holds for all integer m > Kd: 

M (1 -t- A"") < mA"'-^ \\7A\^ , for all \x\^ + \^\^ > 12. 

(c) There exists a constant a2 > 1 such that 

\WA\^<a2A, Vx, ^eM'^. 

(d) For any \x\, |^| > ^ and m > 1, it holds 

A(l + ^"') <mA'"-^ |VA|^ e^6{m^l + 3d). 
Proof. Using Young's inequality we easily obtain 

(4.2) ^i\xf + \e)<A{x,o<\^\' + m\ 

(4.3) ^(ixp + len < ivA(x,op < 3(ixp + 

This yields assertion (c). To show (a), we note from (|4.2p that 
1<A<A", V|a;p > 12. 

Hence, we obtain with 

1 + A'" < 2/1™ < 36 A"'-V^P, 

which is assertion (a). We further note that assertion (b) is a direct consequence 
of (a). Finally, to prove assertion (d), we compute 



A(l-f-yl'") ==mA"-i |VA|^ 



TO - 1 2d 



A \VA\ 

<mA"'-^ \VAf eH(m-l + 3d), for Iccl, 1^1 > -. 

e 
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□ 

Remark 4.3. Note that the constants ai > ^,02 < 36, and K ~ < 144 can 
be chosen independent of m G N and of the spatial dimension d G N. Moreover, 
J-iT = 1 is not possible for d = 1. 

4.2. Semigroup properties on T-Lm- Analogously to p.7p . we now define the 
following closed subset of T-Lm'- 

T~Lm {W e T-Lm ■ W _L Vm} , TO € N, 

which is again characterized by the zero-mass condition 



(4.4) {w,iym)n^= w{x,Odxd^^0. 



Thus we have Ti^ C Tim Vm e N. As before, we define L on Tim via L = L\„^ 
which yields a closed, densely defined operator on Tim for each to e N. In addition, 
we also have the following result. 

Lemma 4.4. For each to G N, the operator L generates a Co-semigroup of bounded 
operators on T-Lm, satisfying 

(4.5) \\Ut\\mn^)<^^"''' 



Proof. We compute 



Rc{Lw,w)n„, - // -dxd^ 



2d Vri 



^11 '"''^ (^AA"-toA™-i |vA|V^j dxd^, 



by taking into account (|4.ip and the fact that i',^^ = 1 + A™{x, £,). Using assertion 
(c) of Lemma 14.21 we can estimate 

A(A™) =toA™-2 ((„^_ 1) \s/A\'^ + 2dA) <m{{m~ 1)02 + 2d) A"'-\ 

Moreover, since, for all x,(, E M'': A'"~^{x,^) < 1 + A™(a;,^), we consequently 
obtain 

A(A") - TO A™-! |VA|^ + — < A(A™) + — < /3,„(1 + A"), 

where 

l3m 2d + to(^(to - 1)02 + 2d^ . 

In summary, this yields 

Re {Lw,w)^^ < Pm 
Thus, for the unperturbed evolution equation dtw = Lw we infer 

^Ihll^^ ==2Re {Lw,w)^^<2Pm\\w\\l^^, 

and the assertion follows. □ 

Remark 4.5. Note that /?„ cannot be negative in Lemma 14.41 since Utifi) — fJ-- 
However, using some refined estimates below, we shall find (see Proposition 14. 8p 
that the restricted semigroup C/t|^± is exponentially decaying, provided to S N is 
sufficiently large. To this end, we note that the two (non-orthogonal) subspaces 
= span[/i] Tim invariant under L and under Ut (Vto G N) due to mass 
conservation p.8p and (|4.4p . 
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As a final preparatory step, we shall need the following decomposition result for 
L, where we denote 

nl, {w e n.n ■■ e -H™}. 

Proposition 4.6. Let m > Kd be some fixed integer, and K was defined in Lemma 
\4-'^ Then there exists an < e < 1 such that the operator L can be split into 
L = L\ + L2, with ifjil defined in HA.1\) and HA . ^|) and satisfying: 

(1) LI : Hm %m is a closed and unbounded operator, while L\ : Hl^ — > H 
and L\ : Hm ^ are bounded operators. 

(2) (LI — z) : % ^ % and (LI — z) : Jim — > 'Hm oltc closed, unbounded and 
invertible operators for every z e 17 {z e C : Rcz > — A,„}, where 
Am > is a positive constant defined in (jA.5|) . 

(3) The operator 

^1(^2 ~ ^) ■ "H-rn T-i (Z "Hm 

is bounded for any z G f2. 

The proof is lengthy and rather technical and therefore deferred to Appendix A. 

Remark 4.7. Note that in ProDOsition l4.6[ e has to be chosen positive, in order to 
ensure assertion ([1]). In fact, while L| continues to be coercive also for e = 0, the 
operators Lf : Hm ^ and Lf : Hm. H^^ become unbounded as e — 7> 0. The 
fact that Lf is bounded for e > is essential, in order to obtain the decay estimate 
(|46)) . cf. the proof of Proposition l48l 

Indeed, introducing the decomposition L = Lf + L| is one of the key ideas in 
mil^n] in order to lift estimates for the resolvent R{z) = {L — z)~^ onto the larger 
space Hm- The general decomposition procedure introduced in [24] applies to the 
WFP equation and provides the following exponential decay of Ut on Hm- 

Proposition 4.8. Let a > be the spectral gap of L^ in H, and let wq G Hm with 
JJ]gi2d Wq dxd^ = 1. Then, for every integer m > Kd, it holds 

(4-6) \\Ut{wn - < SmC^^-^Wwo - 

for anyjm € (0,7,„), where 

(4.7) jm ■= min{A„ ; cr} > 0, 

and dm — Sm{jm) > I is given in (|4.15|) . Furthermore, we have for the resolvent 
set 

(4.8) emnj ^^i-={zeC:Rez> -7,„, z ^ 0}. 

Proof. For the sake of completeness we briefly present the proof which follows the 
ones of Theorem 2.1, Theorem 3.1, and Theorem 4.1 in [53]. The spirit of the proof 
is the following: By using the operator factorization from Proposition 14.61 we shall 
infer an estimate for the resolvent on Hm- Restricting the resolvent to Hm removes 
its singularity at z = and consequently yields a uniform estimate on the complex 
half plane {z e C : Rez > —7m}. The Gearhart-Priiss Theorem [53JI32] then yields 
the exponential decay of Ut on Hm- The proof now follows in several steps: 
Step 1: Following [2H] we define on Hm the operator 

(4.9) i?(z):=(L^-z)-i-((L-z)|^)"^Lf(L|-z)-\ zeVli. 

Theorem 2.1 and Remark 2.2 in [24] implies that R{z) is the inverse operator of 
(L — z) in Hm for any z G fii and therefore statement (j4.8p holds. 

Step 2: The next step is devoted to obtaining uniform estimates for i?(z) = 
^(L — z)|^ ^ , for z G C on some appropriately defined half planes, cf. [24l 
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Theorem 3.1 (4)]. To this end, we shall first prove the following bound for the 
resolvent on 1-L: 

(4.10) sup||(i- (a + js))-i||.^(^) = A'o < oo, Vae(-cr,0). 

s 

Indeed, the constant Kq = Kq (cr, a) can be explicitly obtained by considering the 
resolvent equation for Rez > — fi and z 7^ 0: 

{L - z)f = g on U, 

Using the orthogonal decomposition / = + ci/i, g = + C2/i, having in mind 
that L maps Ti.^ into "H^, we infer 

/ =KL-z)\^g , ci = . 



In view of p.ll[) this yields 

(Rez + a)^^^^ ' |z|2 {{Rez + a)^' |z|2 



Hence, 



\{L-z) ^WsecH) < max | ^ i |~| | for Rc z > -cr, z 7^ 0; 



and Kn < max \ —J— : rr > . 

" — 1^ a+cr ' |a| J 



-1 



From this bound on {L — z)|^ we can deduce a bound on [L — z)|^ for 

z € r^i = {z e C : Rez > -7™, z ^ 0}. Using gH), (|X4| . and if G ^in}^ U), 
we infer 

II (L - z)^i||.^(«„,) 

< ll(-^2 ^ zy^\\m'H^^-Hl,J (1 + IK-^ - ^)""^ll,'^(«) l|-^lll,«(«i,^-H)) 

= :t9™(z). 

Next we consider the resolvent of L\^^. First we note that both subspaces of 
'Hm ~ span[/^] 'K^ are invariant for (L — z)^^ (cf. Remark 14.51) . Hence (L — z)\^ 
and (L — z)|^^ coincide on 'H.'^^. Since z = is an isolated and non-degenerate 
eigenvalue of , we conclude (t(L|^^) = a{L\^ ) \ {0}. Since L generates a 
Co-semigroup on it is closed and its resolvent is analytic on 

e{^\u^ 3 {z G C : Rcz > -7m}. 

For any fixed 7™ G (0,7m) we henceforth conclude from (j4.1ip that the resolvent 
of L|.^^ is uniformly bounded (on a whole right half space) 

(4.12) ||(L- z)-i|[^^^_,^ < A/ (7™) < 00, Rcz > -7™ . 

Note, however, that the constant Af(7„i) is not known explicitly. 

Step 3: Next we shall show that this resolvent estimate yields an exponential 
decay estimate for the semigroup Ut on H^. In order to do so, we will apply the 
Gearhart-Priiss-Theorem to the rescaled semigroup e'''"'*J7t, cf. Theorem V.1.11 in 
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[T^ (see also [211 Theorem 3.1] and [H]). This is possible in view of the uniform 
bound (|4.12p and yields the following estimate for Ut- 

(4.13) \\u,U^^.j<^-^±^^e-~^-\ t>0, 

where w > /?,„ + 7,„ + 1, Cl < 7r(w - /?,„ - 7„i)"\ and 

(4.14) A/ :-sup||(L- (-7™ + is))-i|U(«^) < i?,„(-7m), 

where the second inequality follows directly from (j4.11[) and the definition of M. 
Note that ()4.14p asserts a bound on the resolvent along the (fixed) line 

{z e C : z ^ -7m + is}, 

in contrast to (j4.12p . Keeping this in mind, we conclude that the estimates in the 
proof of Theorem V.1.11 in in fact only depend on the resolvent evaluated at 
z ~ —7m + is. Interpolating (|4.13p with (j4.5p . we consequently conclude 

II C/t II, < <5me~'''"*, 

where (after optimizing in a; > f3„i + "fm + 1) 

(4.15) Sra := max{|i?m(-7„02 ; e^-+^™ } > 1. 
This finishes the proof. 

□ 

Proposition 14.81 implies that the operator L is invertible in the space ; foi" 
m > Kd. More precisely, invoking classical arguments (cf. Theorem 1.5.3 in |30j). 
we infer 

(4.16) ym>Kd: \\L-' \\^^^^^ < ^ ^: ^. 

Note that the constant am > depends on m and thus on d and A. The reason why 
wc obtain exponential decay of Ut on "Hm with a rate 7m < 7m can be understood 
from the fact that in order to apply the Gearhart-Priiss Theorem one needs to 
guarantee a uniform resolvent estimates on some complex half plane, which is not 
sharp, in contrast to e.g. the Hillc-Yoshida theorem (used on H). 

5. BOUNDEDNESS OF THE PERTURBATION 

In this section we shall prove the boundedness of the operator & [Vo] in Tim which 
is the key technical result for our perturbation analysis. We recall that our potential 
V from (j2.1[) consists of the harmonic potential plus the perturbation XVq . Hence, 
we shall now consider — A0[Vb] as a perturbation of L. 

Proposition 5.1. Let m G N and the potential Vq E C™(M''). Then the operator 
Q[Vo] maps H,n into cind 

l|0[^o]||.®(^i„,) < Tm := Cm max \\diVo\\ ^^^^ , 

where Cm > denotes some positive constant, depending only on m and d. 
Proof. We first note that, in view of (j4.2p . the norm ||/|||^^ and the norm 

ll/ll^„.:= Jljf\"(^^^) (l + (N' + m") dxde 
are equivalent. It is therefore enough to prove that 
(5.1) // \e[Vo]w\' (l + (|a:|2 + |^|2)™) d.TdC<Co||Hlm: 
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for some Cq > 0. In the following we denote by 

(5.2) {J'i-.v^Kx, v) = w{x, 77) := / w{x, ^ e-'^-^d^, 

the partial Fourier transform with respect to the variable e M'' only. Recall from 
(|1.3p that the operator Q[Vo] acts via 

(5.3) Q[Vo]w = -iT-^^{SVoix,ri) ■ T^^^wix^r])). 
Using Plancherel's formula and Holder's inequality, this implies 

\\emw\\^.<2\\Voh^ \\wh2. 

Thus, in order to prove (jS.ip . we only need to estimate 

\e[vo]w{x,o\' {\x\^+\er dxde 

We rewrite this term, using 

j=0 ^ 

in the following form: 

= E (7) '"^''^""''^ \emw{x,o\^ dxdc 

It holds 

where Cnj are some coefficients depending only on n G N''. Therefore 

e[Vo]w{x,of {\xf + mT dxd^ 

= E M I E / 1^1'^""^' mvoMx,o\' dxd^ 

where we denote := ^^"^ • • • ^^"'i. From (|5.3p we see that 

(5.4) lir mVo]w)\\l.^^..) = (2^)-^ (.5^0^^) II'. (R..) , 

with w{x,r]) defined by (j5.2p . We expand the right hand side of this identity by 
using the Leibniz formula (see also (U]), and we apply 

sup \dl (SVo) {x,7j)\ < 2(1-^) sup \dl Voiy)\ , 
(cf. Definition (|1.4p ). Then we can estimate (|5.4p as follows: 

e''mVoMx,0\^ dxdi < C(n) niax ||a^yo||i^(K.)||r-'^«lli.(K..), 
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where C{n) > depends only on binomial coefRcients. In summary, we obtain 

|e[yoMx,e)P(N2 + iePrdxde 



2d 

rn 



< max lia^FollioofK.) / / + dxd^, 

b\<m J Jm.^'i 

with Cm, Cm > depending only on binomial coefficients. Thus, the assertion is 
proved. □ 

Remark 5.2. The unboundedness of 8[Vo] in H (the exponentially weighted 
Hilbert space) is due to its non-locality, which can be seen from the following 
reformulation of p.3p (cf. §3 of jj): 

{Q[V]f){x,0 = -2n-''W{x,0 H fi^^O, 

with 

W{x,k) := Im [e^'^-'=l>(2fc)]. 

To illustrate the situation let us take Vo{x) = sm{x ■ kg) for some ko E \ {0}. 
This implies 

{&[V]f){x, = 2^ cos(2a; • fco) [f{x, f - fco) - /(x, ^ + fco)]. 

The problem is that such a shift operator cannot be bounded in an L^-space with 
an inverse Gaussian weight, as the following computation shows: 

6. Proof of the main theorems 

The results of the preceding sections allow us to give the proofs of Theorem [T] 
and Theorem [21 

Proof of Theorem\^ We start with Assertion (i): Let m be the integer fixed in the 
assertion. Any solution Woo G of (j2.2p that is subject to the normalization 
// Woo dxd^ = 1, satisfies the unique decomposition Woo = IJ. + w^, with w^, £ T-Lm, 
i.e. //u;*d.Td^ = 0. Therefore, we consider the following fixed point iteration for 

T : Urn T-Lm, It'n-l ^ T{Wn-l) = Wn, 

where Wn e Hm solves 

LWn = A e [Vo] {Wn-l + H) ■ 

To be able to apply Banach's fixed point theorem, we have to prove that the map- 
ping T is a contraction on Hm- To this end we write, for any Wn-i,Wn-i G 

II - Wn\\-H±^ = II Ai"^ e [Vo] {Wn-1 ~ Wn-l) ||^_l 

and estimate 

\\Wn - Wn\\f^± < |A| II L^^ \\s8{-H^) l|0[Vb](w„_i - Wn-l) \\fi± ■ 

From (|4.16p and Proposition 15. II we obtain 

^ TmlAI _ 
\-U± < ||W„_1 - 'Wn-l\\-f^± , 
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since the potential Vq satisfies (|2.9p . Since |A| < (Jm/^m, there exists a unique 
fixed point = T{w^) € H^j. Thus, the unique (stationary) solution of (j2.2[) is 
obtained as Woo = /i + G Hm. Note, however, that /x / in the sense of Hm. 

The obtained solution Woo is real valued, since T maps real valued functions to 
real valued functions. Moreover, Woo G Hm satisfies (|2.2p . at least in the distribu- 
tional sense. Furthermore, 0[Vb]iUoo € Tim and Lwoo G 7J~^(R^'') and thus (j2.2p 
also holds in (R^'') . To explore - a posteriori - the regularity of Woo , we rewrite 
(j2.2p in the following weak form 

// (V^eWoo • Vj;V3 + V^iuoo ■ V5(/7 + 'u;ooV')dxdC = H-^{F{woo),^)m, 

J JS?d 

for any G H^i^'^), where 

F{woo) ■= Woo - diva;(^Woo) + div^{xWoo + S^u^oo) - A6[V'o]woo. 
Clearly F(woo) G i?-^(K2d) and thus Woo G H'^{R^'^) follows. Moreover, since 
F{woo) G we also have Woo G ^[^^^(I*^'')- 

For the proof of Assertion (ii), we first note that &[Vo] is a bounded perturbation 
of L on Hm and thus (jl.ip admits a unique mild solution w G C([0, oo), "Hm)- 
Since 0[Vb] niaps into "H^, we also know that along this solution the mass is 
conserved, i.e. // w{t) dxd£, = 1, for all t>0. 

Next, consider the new unknown g{t) := w(t) — Woo with go — wq — Woo- Due to 
mass conservation g{t) G Hm for all i > 0, and we also have 

dtg = Lg-XQ[Vo]g, 

since Woo is a stationary solution of (|l.ip . Taking into account that the semigroup 
Ut associated with L in the space Hm satisfies (|4.6p . it holds: 

Umn^. < Sme-'^-'WgoWn^ + Sm\X\ f e-^^-^'-'^mVoMs^n^ ds 

Jo 

< Sme-'f-'hoWH^, + SmlMTm f C"^™ t*"^) 1 1 <7(s) 1 1 ds, 

Jo 

for any 7,„ G iSm\X\TrmJm,) and with Tm defined in Proposition 15.11 Gronwall's 
lemma then implies 

||5(<)ll«,„ <5rnC-*(^'"-*'"l^|I^'")|1.90|l«„. 



It remains to prove Assertion (iii): As before we write Wac = M + w^, where 
G Hm solves 

{L-xe[Vo])w, = xe[Vo]fi. 

Since ||0[Vb]|l ^(-^^j < r„i, we obtain || A9[Vb]/-t||.^^ < |A|r,„j|/ij|-H,„ . Next, we 
consider L on H:^^. From the proof of Proposition 14.81 we conclude for its resolvent 
set: 

Q [L\'H± ) 3 {z G C : Rez > -7™} , 



and thus 

g (^{L - Xe[Vo])\^^^ D{zgC:Rcz> |A|r„,-7„}. 
Since |A|rm — 7m < (see (|2.10p ) we have 

{L - Ae[yo])"' = L-^ (Id - Xe[Vn]L-^y' (on Hi). 



STATIONARY STATES AND LARGE TIME BEHAVIOR OF WFP 



19 



Using (|4.16p and |A|r„i < cr,„ (see (|2.10|) ) we conclude 



\\{L-XQ[Vo]) 
Thus, by writing 



< 



1 - |A|r„i^ cr„,. - |A|r„ 



w,^({L- xe[Vo]) 



we infer 



I'M,, 



Wo 



Ml 



Hi. 



< 



{xe[Vo]fi), 

\X\T„M\n^ 

Cm — A r,„. 



and the assertion is proved. 



□ 



Remark 6.1. Due to the mass normahzation J J Woo da;d^ ^ // pdxd^ = 1, the 
fixed point must take both positive and negative values. Thus, Woo ~ p. + w^, 
may, in general, also take negative values. 

Proof of Theorem We start with assertion (ii) , which follows from the fact that 

II P 11.^, = II Pi; •) IIl. = i'^n)'/' II •) IIl^ < (2^)'/' II ■)\\n„.^ Vp e 
Thus, we infer 

p{t) Poo in , 
with the exponential rate obtained from Theorem 1 (ii). 

To prove assertion (i) we consider the transient equation (jl.ip as an auxiliary 
problem: Choose any po G ST'^ such that tr po = 1 and the corresponding wo € "Hm . 
Due to the results on the linear Cauchy problem given in [5] we know that p.ip 
gives rise to a unique mild solution p £ C([0,oo); satisfying ti p{t) = 1, for 

all t > 0. Hence, the trajectory {p{t), t > 0} is bounded in Since 5^ has a 
predual, i.e. the compact operators on L^(R'*), the Banach-Alaoglu Theorem then 
asserts the existence of a sequence {tnlnGN C R+ with i„ — > oo, such that 



P(tr. 



in weak-* 



for some limiting p G The already obtained .^-convergence of p{t) towards 
Poo G ■% implies poo = P & ^i- And the uniqueness of the steady state yields the 
convergences of the whole t-dependent function p{t) — >■ p in ^ weak-*. Finally, 
we also conclude positivity of the operator p^o by the ,52-convergence and the fact 
that we already know from [S]: p{t) > 0, for all t > 0. 

It remains to prove tr p^o = 1. To this end, we recall that for any p £ the 
corresponding kernel 

d(x,r]) := p(^x+^ , x-^^ 
satisfies € C{R^, L\{mi)), see [J, and it also holds 



.1) 



tr p 



■d{x,0)dx. 



Further, note that i!){x, r/) = {T^^jjw){x, rf) = w{x, r/), by p.2p . On the other hand, 
for any w £ Hm we know that w £ C(M^, i^(M!^)), due to the polynomial L^-weight 
in X S K'* and a simple Sobolev imbedding w.r.t. the variable 77 G M'' (for both 
embeddings we used m > ^). Hence the normalization condition jjwao dxd^ = 1 
implies trpoo = 1; via (j6.ip . and assertion (ii) is proved. 

Finally, we prove claim (iii) by first noting that the ^^-convergence of p{t) implies 
convergence in the strong operator topology. Thus, having in mind that ||/3(t)||^^ = 
IImoII^^j = 1, we infer from Griimm's theorem (Th. 2.19 in [35]) that p{t) also 
converges in the ^-norm towards p^o- This concludes the proof of Theorem[21 □ 
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Appendix A. Proof of Proposition 
The proof will be divided into several steps: 

Step 1: Let x e C^(]R^'') be such that x = 1 on Si(0), with supp(x) C ^2(0), 
IIVxIIl- < 72, and let Xe{y) ■= x{ey). for any y = {x,0 & < e < 1. We 
define 

(A.l) Lfw := (dw - j/^ Vt-"^ • Vw) Xe, 

as well as 

(A.2) L^w := I'm div Vw) +S/wS/A + dw + div{w F) 

+ {dw- iy„i ^v~l ■ Vw) (1 - Xe) • 

It is easily seen that L\ indeed satisfies property (f ). 

Step 2: In order to prove properties (2) and (3), we need to show that 

{LI — z) ^ : Hm — > V-ln^ 
is bounded for z £ C C. To this end, it suffices to show that (Lf ~ ^) satisfies 

Wz e n -. 

(A.3) -Re {{LI ~ z)w, w)n^ > c\\wf^, , 



with some c ~ c(Rez) > 0. Indeed, suppose (jA.Sp holds for all (complex valued) 
w G -.^ {w e n,n ■■ Vw, (V^ + F) ■ Vw,Aw G nm}- Then, L| - z is densely 
defined on i^(L|) := Ttf^ C Hm and dissipative. Hence, it has a maximal dissipative 
(and thus surjective on Hm) extension, which we shall consider in the sequel. From 
(IA.3I) we conclude 



c\Ml^i^ <-Re{{Ll-z)w,w)-H^ < \\w\\ni,J\{Ll - z)w\\n^, c> 0, 
which implies 

(A.4) ||(L| - z)-hd\\ni^ < I fo'- all ^ E Hm- 

For future reference we note that this bound will turn out to be uniform on the 
lines z = a + is, s £ M. (with fixed a > — A„i), since c = c(Re z). 

Step 3: Now, we have to prove (|A.3|) . To this end, we decompose 
Re {Llw,w)^^=- JJ (l + A")|Vwp dxd^ 

|^|2^m-l |y^|2 ^^^^ 



m 
2" 



\w\^ (1 + A") (1-Xs) dxd^ 
|wpV(l-Xc)-V(l + ^™) dxd^ 
(1-Xe) A(l + A™) dxd^ 



1 

+ 2 
1 

+ 2 

h + l2 + l3 + h + h 
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Us: 



ing Lemma 1321 (b) , we readily obtain for e < 



m 



In order to treat the term I4, we note that 

iVxel = £|Vx|, |Vx| < V2, supp {Vxe} C {y e K^"^ I J < |y| < J}. 
Therefore, by (03), 

< yiSe for y = (a;,^) e supp {Vxe}. 

With £ < 1 this aUows us to estimate 

|V(l-x.)-V(l + A™)| <£mA™-i |Vx| |VA| 

which imphes 

I/4I < 3me // IwfA""-^ \VAf dxd^. 



The term I5 can be estimated using Lemma [421 (d): 
h^l II (1-x.) A(l + ^™)dxde 



2 _ 

< Sme^ (m - 1 + ad) // |wp (1 - Xe) ^""M VA|^ da; d^ 



In summary, we obtain 

i^^ + l^^l^ V8 



^/2 + |/4|<-fv-3™e) // Iz^r^^^MV^r dxdC, 



and 



^/2+/3 + /5 <-(v-3me2(TO-l + 3d)) // \w\^ A""-^ \V A\' dx dt 
4 V 8 / J jR2d 

Now choosing e < min 12^7^ j"j '^6 can estimate (using m> d) 
™ ( ^ - 3e ) > 0, ™ f ^ - 3£^(7n - 1 + 3d) ) > ^. 
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Therefore 



-Rc (L>, > // (1 + A'")|V?i;pda;d^ 

J JM2d 

+ (9-2) ffjw\'A^-'\VA\'dxd^, 



which we estimate further using Lemma 14.21 (b) : 

-Rc {LIw,w)h^> II (l + A'")|Vu;p dxdC 

J JK2<i 



H'(l + A'")(l-Xi/vT2) da:de 
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This establishes the desired estimate outside of {x, ^) G B^^^^{0). In order to take 
into account the contribution near \y\ = 0, we consider 

i 

Applying Sobolcv's inequality we obtain for any d > 1: 

(l + A"-) \Vw\' dxdC > [[ \\/wf dxdC > 



with q = 21=2 ■ This can be estimated further via 



(J2 

Cd l|w'lli,(R2.) > Cj ||w||i,(s^^^(o)) > 111 + A"^\\L^^B,^^,m "'""«™(-Bi/vT2(0))' 

where C'l depends only on the Sobolev constant Cd and on the measure of the ball 
^i/^/ni^)- Finally, in order to deal with d = 1 we apply Cauchy-Schwarz to obtain 

\Vw\dxd^= [[ |Vw| + (l + A™)^^/^da;d^ 

- (// '^""'^ (1 + ^"')^^'^^) ^ (//, (l + ^™)"'da;de^ ^ . 

By assumption we have m > 1 (cf. Remark l4.3p . Hence, the second factor on the 
r.h.s. is a finite constant, denoted by CA,m- Applying again Sobolcv's inequality 
(for d = 1), we obtain 

~ C 

||V«;|Ui(K2) > Ci \\w\\l2(r2^ > — -^yi ||w||«„(s^^^(o)). 

By combining all the above estimates we infer for all (complex valued) w G H^'- 

—Re {L2W,w)^^ > A,n {w,w)^i , 

where Am > is given by 
(A.5) 

. r Cl Cl 1 

In summary, inequality (|A.3p holds on "H^n for any Rez > —A„i. And we easily 
see that c{z) = min{A,„ ; + Rcz}. Hence, the operator (i^ ~ z) is invertible in 
T-Lm- And we have to choose e = e{m) > such that 



(A. 6) e < min 



1 1 
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